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Traditionally the effects of MHD instabilities and micro-instabilities on plasma confinement 
are investigated separately. However, these two instabilities often occur simultaneously, with 
^ . the overlap of the dynamics on a broad range of spatial scales. It is thus vital to incorporate 
^ I these instabilities consistently by a proper multi-scale modeling. Furthermore, there has been 
an overwhelming evidence that the overall transport of heat and particles is significantly influ- 
enced by intermittency (or bursty events) caused by coherent structures. A crucial question in 
plasma confinement is thus the prediction of the probability distribution functions (PDFs) of the 
PL|. transport due to these structures and of their formation. 

In this paper, we investigate intermittent transport in a multi-scale model by consistently in- 
■ corporating both tearing instabilities and micro-instability due to pressure gradient [IJ. We first 
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present an exact nonlinear solution in the form of a coherent structure (modon or bipolar vortex 
soliton). We then present a first analytical result of intermittency in our multi scale model by 
utilizing a novel non-perturbative method. Specifically, we compute the PDF tails of momentum 
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flux and heat flux, by assuming that a short-lived modon is a coherent structure responsible for 



^ . bursty and intermittent events, contributing to the PDF tails. 
> 

' The governing equations for modeling the MHD and drift wave turbulence [IJ are as follows, 
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^ + [(^:P] = ~v,(^{l-K2)^ + K2^^+C^[xir,Vlxir]+XiVlp+f2. (2) 
^ + [^-p,xjf] = -v^^ + nVlxj/ + h. (3) 



■ The model is a reduced MHD model that includes tearing modes and interchange modes and 
consists of a set of coupled equations for the electrostatic potential ^, the pressure p and the 
magnetic flux Xj/. Eq. (1) is the combined electron and ion momentum balance equations where v 
is the viscosity and Eq. (2) represents the energy conservation (neglecting parallel ion dynamics) 
with Xi as the diffusivity. Eq. (3) is the Ohm's law with the resistivity rj . The magnetic field is 
assumed to be dominated by a strong component in the z-direction. Here the constants are = 
^ ~ B^/(2^i)) ' pressure gradient, L± is related to island width, f2 = ^ is the ion 

cyclotron frequency and Ta is the Alfven time. Eqs. (l)-(3) are normalized by using the Alfven 



speed va = L^/ta- The interchange mode is controlled by curvature and pressure gradients 
through the parameters K"i = 2Q/Ta^ and K2 = {Rq is the major radius), respectively. The 
tearing mode is determined by the parameter C = jjj-Q.fmi. The dynamics of magnetic islands 
are strongly influenced by the magnitude of C. The statistics of the forcings ? = { 1 , 2, 3}) is 
assumed to be Gaussian with a short correlation time modeled by the delta function as 

{fi{x,t)fj{x\t')) = d,jd{t-t')Ki{x-x'), (4) 

and (//) = 0. The angular brackets denote the average over the statistics of the forcing fi. The 
delta correlation in time is chosen for the sake of simplicity of the analysis. 

For computing PDF tails the basic idea is to associate the bursty event causing intermittency 
with the creation of a coherent structure (e.g. blobs, streamers etc.). The creation process is 
identified with the instanton which is localized in time, existing during the formation of the 
coherent structure. This idea is embedded in the instanton method, which is a non-perturbative 
way of calculating the PDF tails. We calculate the PDF tails of local flux M{^{x = xq)) that 
is the second moment of ^ (e.g. momentum flux or heat flux). The PDF tails are expressed in 
terms of a path-integral using the Gaussian statistics of the forcing [4J . The optimum path is then 
associated with the creation of a short lived coherent structure (among all possible paths) and 
the action is evaluated using the saddle-point method on the effective action. The probability 
density function of the flux R can be defined as 

P{R) = {5{Mi^{x = xo))-R)) = JdXe'^\, (5) 

where 

I^ = {exp{-iXM{(j){x = xo)))) = j &<^&^&p^p^\if^^fe^^K (6) 

We assume that the instanton solution is of the form = F{t)^{x,y), p = G{t)p{x,y) and 
W = J{^)x{^:y) where j = {cp.p.x} are the modon solution satisfying V^Xj = OCijXj + a2jX. 
Furthermore, we assume a linear relation between and p. To simplify the computation of 
the path-integral the saddle-point method will be used where the variational first derivatives in 
((^, ^, p, p, i//, \fr) of the action will be put to zero. The saddle-point solution of the dynamical 
variable ^(xj) of the form ^{xj) = F{t)(p is called an instanton if F{t) = at ? = —0° and 
F{t) 7^ at ? = 0. Note that, the function (p{x) here represents the spatial form of the coherent 
structure. Two limiting cases will be studied in detail; (1) iso-instanton case {F = G = J); and 
(2) the electrostatic limit (F = G and 7 = 0). In the iso-instanton case (1) we use forcings in 
Eqs. (l)-(3) whereas in the electrostatic case we use forcings in Eqs. (1) and (2). Note that the 
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Figure 1: (Color online). The ratio of the coefficients ^2 (electromagnetic) and ^2 (electrostatic) 
as a function of the parameter C, for the other parameters K"i = 0.01 and K2 = (blue line, 
dashed line) and k"i =0.1 and K2 = (red line, solid line) are shown. 

conjugate variables that mediate between the physical variables and the observable flux arise 
due to the uncertainty in the physical variables of the stochastic forcing. 

In the iso-instanton case the PDF tails of momentum flux and heat flux are found to be 
P{R) ~ \fRe^^'^^^^ and for particle density P{J^ ~ R^e^^^^^ . Here the saddle point method 
with movable maxima [I6JI has been used to compute the integrals in Eqs (5) and (6). Although 
the power laws multiplying the exponentials have not been included in previous studies, the 
exponential scaling were obtained previously. The tails of PDF of momentum flux and heat flux 
are found to be stretched exponentials which are broader than a Gaussian. This suggests that 
rare events of large amplitude due to coherent structure are crucial in transport (similarly to 
what was found in the previous studies f2\, |[5l), offering a novel explanation for exponential 
PDF tails of momentum flux found in recent experiments at CSDX at UCSD dSJ. These and our 
previous results highlight the key role of structures on intermittent turbulent transport. 

In a similar way the PDF tails in the electrostatic case is computed and the PDF tails of 
momentum flux and heat flux are found to be P{R) ~ ^/Re^^^'^^^^ and for particle density flux 
P{R) respectively. 

In particular, we will now examine the crucial dependence of the overall amplitude of the 
PDFs on the relevant physical parameters. We show the effect of magnetic fluctuations on mo- 
mentum transport by an explicit comparison of the coefficients in the iso-instanton {^2) case 
with the electrostatic case (^2) in Figure 1. It is shown that the values of coefficient in the elec- 



tromagnetic case is larger than those found in electrostatic case in a scan of the parameter C for 
two sets of parameters fCi = 0.01 and ^2 = (blue line, dashed line) and fc = 0.1 and K2 = 
(red line, solid line) (the parameters are similar to those used in Ref [IJ). Note that by using 
K"2 = the linear diamagnetic effects are suppressed and = and that viscosity (v), diffu- 
sivity {x±) ^rid resistivity (rj) are taken small to be ~ 10^^. This suggests that the PDF tails in 
the electromagnetic case is smaller and thus the effects of rare events are smaller in comparison 
to the electrostatic case. Note that the parameter C determines the temperature response of the 
tearing modes. Note also that the modon structure is determined by matching the inner and outer 
solution in the same manner as the drift wave vortex solution in Ref. [171. 

Interestingly, there exists a modon solution where the Maxwell stress is canceled by the 
Reynolds stress, resulting in a ratio ^2/ <^2 ~ 1 for all values of C, Ki and K2. This illustrates the 
exact cancellation of Maxwell stress by Reynolds stress due to electromagnetic fluctuations. In 
the present study the effects of pressure gradient are small since the parameter K2 = 0. 

We note that treating first the MHD and micro-turbulence as a fully multi-scale model as in 
Ref. JHl and then computing the PDF tails could possibly lead to other exponential dependencies 
of the PDF tails. This will be investigated in future publications. 
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